Dynamic response of strongly correlated Fermi gases in the quantum virial expansion 



O 

(N 

wo: 
< 

(N 



C/2 

c ■ 

c3 ■ 

+-* ■ 



I 

c 

o 
o 



> 
o 
\o 
o> 
cn 

od 
o 
o 



X 



Hui Hu x ' 2 H and Xia-Ji Liu 1 , and Peter D. Drummond 1 

1 ARC Centre of Excellence for Quantum-Atom Optics, Centre for Atom Optics and 
Ultrafast Spectroscopy, Swinburne University of Technology, Melbourne 3122, Australia, 
2 Department of Physics, Renmin University of China, Beijing 100872, China 

(Dated: August 25, 2010) 

By developing a quantum virial expansion theory, we quantitatively calculate the dynamic density 
response function of a trapped strongly interacting Fermi gas at high temperatures near unitarity. 
A clear transition from atomic to molecular responses is identified in the spectra when crossing from 
the BCS to BEC regimes, in qualitative agreement with recent Bragg spectroscopy observations. 
Our virial expansion method provides a promising way to solve the challenging strong-coupling 
problems and is applicable to other dynamical properties of strongly correlated Fermi gases. 
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INTRODUCTION 

The dynamic structure factor (DSF) plays a funda- 
mental role in understanding quantum many-body sys- 
tems [lj: it gives the response of the system to an ex- 
citation process that couples to density. Experimental 
advances in Bragg spectroscopy have now measured the 
DSF for strongly interacting ultra-cold fermions. At a 
high momentum transfer, these experiments show evi- 
dence for a clear transition from atomic to molecular re- 
sponse when traversing from the BCS to BEC regimesQ. 
However, while the DSF for weakly interactingbosons is 
well-studied theoretically and experimentally [3j, the ex- 
pected dynamic structure factor of strongly interacting 
fermions remains unknown, except for a dynamical mean- 
field approximation 0, [B| . This is due to the notorious 
absence of a small parameter for strongly interacting par- 
ticles. The Bragg experiments create a theoretical chal- 
lenge: how can one reliably calculate the dynamic struc- 
ture factor in this strongly interacting regime? Similar 
challenges also arise when one tries to understand some 
recent measurements on radiofrequency spectroscopy of 
strong interacting fermions (B-Q. 

In this paper, we present a systematic study of the dy- 
namic structure factor of a normal, trapped and strongly 
interacting Fermi gas at high temperatures. This is 
achieved by developing a quantum virial expansion for 
dynamical properties of many-body systems. Our ex- 
pansion is applicable to arbitrary interaction strengths 
and has a controllable small parameter. The fugacity 
z = exp(fi/ksT) <C 1 is small since the chemical poten- 
tial \x tends to — oo at large temperatures T. We note 
that dynamical studies using virial expansions have been 
carried out for classical, untrapped gases However, 
previous virial expansion studies of quantum gases were 
restricted to static properties only [h1 [lH . By compar- 
ing the virial expansion prediction 1 121 with the experi- 
mentally measured equation of state [13], we find a wide 
applicability of the expansion for a trapped Fermi gas: it 
is valid down to temperatures as low as OATp (see Ref. 



[12| for details). Here, Tp is the Fermi temperature of a 
trapped ideal, non-interacting Fermi gas. 

Our main results may be summarized as follows. We 
find a smooth transition in the dynamic structure fac- 
tor, from an atomic response to a molecular response as 
the interaction strength increases (Fig. 1). This feature 
agrees reasonably well with recent experimental measure- 
ments although the latter was carried out at lower tem- 
peratures. We show that the spin-antiparallel dynamic 
structure factor provides the most sensitive probe for 
molecule formation (Figs. 2 and 3). The static structure 
factor is also obtained as a subset of our results (Fig. 4) . 
These predictions are readily testable experimentally. 



QUANTUM VIRIAL EXPANSION OF DYNAMIC 
STRUCTURE FACTOR 



We start by constructing the virial expansion for the 
dynamic structure factor 5(q, uS), the Fourier transform 
of the density-density correlation functions at two differ- 
ent space-time points. Consider a harmonically trapped 
atomic Fermi gas with an equal number of atoms (N/2) 
in two hyperfine states (referred to as spin-up, a =t, 
and spin-down, a —i), where Sff(q L ,u) — Sn(q, lu) 
and Sfj,(q, w) = S^ffq, u)). To calculate the total dy- 
namic structure factor 5(q, lS) = 2[S^(q, £j) + 5-|-j,(q, uS)\, 
it is convenient to work with the dynamic susceptibil- 
ity Q, x^'( r , r '; r ) = -(T r n„(T,T)n ( r>(T',0)), where 
ha (r, t) is the density (fluctuation) operator in spin 
channel cr, and r is an imaginary time in the interval 
< r < /3 = l/k B T. 

The DSF is then obtained jl] from the Fourier com- 
ponents Xcra' (r, r'; iuj n ) at discrete Matsubara imaginary 
frequencies iu n — i2nnkBT (n = 0, ±1, ...), via analytic 
continuation and the fluctuation-dissipation theorem: 



S aa > (r, r';ui) = - 



Im Xacr' (r, r'; iu n -> u + i0 + ) 
tt(1 - e-P") 



A final Fourier transform with respect to the relative spa- 
tial coordinate r — r' leads to S aa i (q, ui). 
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Our quantum virial expansion applies to the dynamic 
susceptibility Xo-o-' (r, r'; r > 0), which is formally ex- 
panded as: 

Tr [e-^ H -^e UT h a (r) e - Hr n a , (r')] 
Xcrcr ' = ^ le -p(U-nM) ' ' > 

At high temperatures, Taylor-expanding in terms of 
the powers of small fugacity z = exp(/i/ ksT) <C 1 
leads to Xaa> (r/jr) = (2X1 + z 2 A 2 + ■■■)/{! + 
zQi + z 2 Q 2 + • • • ) = zXt + z 2 (X 2 - X 1 Q 1 ) + • • • , 
where we have introduced the cluster functions X n = 
- ^[e-^e^ifrje-^vtr')] and Q n =Tr n [ e -^], 
with n denoting the number of particles in the clus- 
ter and Tr n denoting the trace over n-particle states of 
proper symmetry. We shall refer to the above expan- 
sion as the virial expansion of dynamic susceptibilities, 
Xaa> (r, r'; r) = zx*<r',i (r, r'; r) + z 2 x<j<j> .2 0, r'; t) + • • • , 
where, 

Xaa',i (r,r';r) = Xi, 

X CTCT ',2(r,r';r) = JCa-JCiQi, etc. (2) 

Accordingly, we shall write for the dynamic structure fac- 
tors, Sao> (q, Oj) = zSaa',1 (q, w) + z 2 S lTtT >,2 (q, w) -I • 

It is readily seen that a similar virial expansion holds for 
other dynamical properties. As anticipated, the deter- 
mination of the n-th expansion coefficient requires the 
knowledge of all solutions up to n-body, including both 
the eigenvalues and eigenstates. Here we aim to calcu- 
late the leading effect of interactions, which contribute 
to the 2nd-order expansion function. For this purpose, it 
is convenient to define Axaa',2 = {Xaa',2}^ = {X 2 } (/) 
and &.S aa \2 = {Saa' ,2} ■ The notation {} means 
the contribution due to interactions inside the bracketed 
term, so that {X 2 }^ — X^—X^p, where the superscript 
"I" in X 2 ^ denotes quantities for a noninteracting sys- 
tem. We note that the inclusion of the 3rd-order expan- 
sion function is straightforward, though involving more 
numerical effort. 

To solve the two-fermion problem, we adopt a short- 
range 5-wave pseudopotential for interactions between 
two fermions with opposite spins, in accord with the ex- 
perimental situation of broad Feshbach resonances. In an 
isotropic harmonic trap with potential V (r) = muj^r 2 /2, 
the solution is known [11]. Any eigenstate with en- 
ergy Ep = e p i + e P 2 (P = {pl,p2}) can be separated 
into center-of-mass and relative motions, $p(r, r') = 
ip p i([r + r'] /2)ip p2 (r — r'). The center-of-mass wave func- 
tion is not affected by interactions, according to Kohn's 
theorem: it is simply the single-particle wave function 
of a three-dimensional isotropic harmonic oscillator, but 
with mass M = 2m. 

For the relative wave function with a quantum num- 
ber p2 = {n p ,l p }, only the branch with zero rel- 
ative angular momentum (l p — 0) is modified by 



interactions. The relative energy e p2 = {2v Hp + 
3/2)fej is determined by y/2T(-v np ) /T{-v np - 1/2) = 

aho I a, where a% = ^/h/mu!o is the characteris- 
tic oscillator length of the external trap potential, 
and a is the 5-wave scattering length. The rel- 
ative wave function is then given by ip p2 (x) = 

A np exp(-x 2 /4 a L)r( T ^ ? )C/(-^ p ,3/2; a ; 2 /2 a L) I with 
A rip being the normalization factor. Here, T and U are 
the Gamma function and confluent hypergeometric func- 
tion, respectively. Other branches with l p 7^ (together 
with the non-interacting counterpart il> p2 \x) for all l p ) 
are given by the standard single-particle wave function 
of a harmonic oscillator, with a reduced mass to/2. 

With this backgrounds, we turn to consider the 2nd- 
order expansion function for the dynamic susceptibility, 

A W , 2 = - {Tr n [e-^e^ha (r) e-*rv (r')j } W 

The trace is calculated by inserting the iden- 
tity J2q\Q)(Q\ = 1- We find A X ^,2 = 
~ Ep,q { e -^P+r(iM Q ) ( P \n a \ Q) (Q \n a ,\P)} {I) . _ 
The sum is over all the pair states P and Q with 
energies Ep and Eq. Expressing the density opera- 
tor in first quantization: (r) = J^. S (r — r^) and 
(r) = ^2jS(r — Vji), it is straightforward to show 
that, 

Ax„<,,2 = -Y,{e-P B ^ P -E Q]c rQ (rr0 j« ; (3) 

where C^ Q = J dv 2 dr' 2 [$^$ Q ] (r, r 2 ) [$^$ P ] (r', r' 2 ) 

and C^ Q = J dridr 2 [$J,* Q ] (r,r 2 ) [$^$ P ] (n,r'). The 
dynamic structure factor can be obtained by the analytic 
continuation, giving the result that AS a <j',2 (r, r'; w) — 

Ep, q {S (u> + E p - E Q ) e-^C P J (r, r')}^ Apply- 
ing a further Fourier transform with respect to x and in- 
tegrating over R, we obtain the response AS ff<T i 2 (qj w )j 

AS^, 2 = + £ P - E Q) ^F P J (q)} (/) , 

P,Q 

(4) 

where F^j? (q) = / dvdv'e-^^'^C^, (r,r'). 

To proceed, one notices that F P ® (q) can be sepa- 
rated into center-of-mass and relative motion parts. We 
thus introduce 5 (cj + E P - Eq) e~ pEp = J du)'5(u' + 
(pi - (qi)S(io — lu' + e p2 — e g2 )e _/3 ^ pl+ep2 ^ to rewrite 

AS roV = {j dw'W cm (q,w')W r ^' - uj')}^ , where 
Wcm = E P i q J(^' + e pl -e ql )e-^ \f plql \ 2 and W$ = 
E P 2 g 2^ - W + e p2 - e q 2)e~^A^' q2 , with f plql = 
JdRe-^^R^R), All q2 ^ \A p2q2 \ 2 and A% 2 ee 
A 2 p2q2 , and A p2q2 = J dxe-^/^M^M- At high 
temperatures, we may apply a semi-classical Thomas- 
Fermi approximation for the calculation of W cm . In con- 
trast, has to be summed over all the eigenstates 
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since the relative wave functions could be spatially sin- 
gular due to strong interactions. 

After some algebra, we find W cm = 
Byjm/it exp[— j3(u)' — ljr/2) 2 /(2wr)], with a con- 
stant B = {k B Tf/ 2 /{qtfiul) and U) R = hq 2 /(2m) being 
the recoil frequency for atoms, and 



A a j' 2 = (-1)^-^(21 + 1) 



dxx 



''jl (y) 4>n p l p (X) 4>n q l q 0) 



,(5) 



where we specify p2 — {n p l p } and q2 = {n q l q }, and I = 
max{Zp, l q }. Here, ji is the spherical Bessel function and 
tf> is the relative radial wave function that can be obtained 
from the two-fermion solution. We require that either 
l p or l q should be zero (i.e., min{Z p ,Z g } = 0), otherwise 
Apf 2 will be cancelled exactly by the non-interacting 
terms. Inserting the expression for W om and Eq. ([5]) into 
AS„^.2, we finally arrive at (w = u) — ljr/2), 



ASu ff '.2 — B 



P 2q2 



2 ) 



-/3e P 2 A era' 1 
A p2q2 r 



(I) 



(6) 

Eq. © is the main result of this work. Together 
with the non-interacting structure factor at large T, 
S l 3,i = C'BvWNe-«-^) 2 /( 4 ^» and 5^,2 = 

-<5<T<r' J B\/m/(167r)e-' 3(tJ - tjR)2/(2 " n) ( 1 + e _/3a 0, we calcu- 
late directly the interacting structure factor, S aa i (q, uj) = 

zS^ t l + z 2 [AS aa ^ 2 + <?^, )2 ], once the fugacity z is de- 
termined by the virial expansion for equation of states. 



RESULTS AND DISCUSSIONS 

Considerable insight into the dynamic structure factor 
of a strongly correlated Fermi gas can already be seen 
from Eq. ([6]) , in which the spectrum is peaked roughly at 
w .R,m = Wj{/2, the recoil frequency for molecules. There- 
fore, the peak is related to the response of molecules with 
mass M = 2m. Eq. ([6]) shows clearly how the molecular 
response develops with the modified two-fermion energies 
and wave functions as the interaction strength increases. 
In the BCS limit where AS aer ' t 2 is small, the response 
is determined by the non-interacting background that 
peaks at u>r: see, for example, S^, l and S^, 2 . In the 
extreme BEC limit (a — » + ), however, AS acr >.2 domi- 
nates. The sum in ASa-cr',2 is exhausted by the (lowest) 
tightly bound state <j) re i(x) ~ y/2/ 1 ae~ x l a with energy 



-h 2 /(r 



The chemical potential of 



molecules is given by fi m — 2/i + ts ■ Therefore, the dy- 
namic structure factor of fermions takes the form, where 
is the molecular fugacity of 



S. 



BEC 



M 

z m B\l — exp 



/3(u - UJR.mf 



4.uj 



R,m 



(7) 
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Figure 1: (Color online). Evolution of the density struc- 
ture factor with increasing interaction strength 1/kpa at 
T = 0.5T F (a) and T — Tp (b). The inset in (b) shows 
the peak position as a function of 1/fcpffl. 



This peaks at the molecular recoil energy uiR ym . As 
anticipated, Eq. ([7]) is exactly the leading virial ex- 
pansion term in the dynamic structure factor of non- 
interacting molecules (c.f. 5^, ,). It is clear that 
5f-|-(q, to) ~ 5-1-4. (q, w) in the BEC limit, since the spin 
structure in a single molecule can no longer be resolved. 

To understand the intermediate regime, in Fig. 1 we 
report numerical results for the total dynamic structure 
factor as the interaction strength increases from the BCS 
to BEC regimes. In a trapped gas with total number 
of fermions N, we use the zero temperature Thomas- 
Fermi wave vector kp — (24iV) 1 / 6 /a/ lo and temperature 
Tp = (3N) 1 ' 3 fiwo/fcs as characteristic units. In accord 
with the experiment [2], we take a large transferred mo- 
mentum of hq = 2>hkp. A smooth transition from atomic 
to molecular responses is evident as the interaction pa- 
rameter 1/kpa increases, in qualitative agreement with 
the experimental observation (c.f. Fig. 2 in [2] ). As 
shown in the inset, the transition shifts to the BEC side 
with increasing temperature. 

According to Eq. © , the molecular response in A5i--|- 2 
and A 5-1-4^ 2 should have the same order of magnitude. 
However, the response is less obvious in the spin-parallel 
channel because of the non-interacting background in 
5ff(q,w). As a result, 5^(q, oj) ~ z 2 A5^ i2 provides 
an ideal probe for the formation of molecules. Experi- 
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Figure 2: (Color online). S-f-|-(q, ui) (dashed lines) and 
iS-f^q, id) (solid lines) at different interaction strengths and 
at T = 0.57> (left column) and T = T F (right column). 
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Figure 3: (Color online). Temperature dependence of 
Sfi(q,oj) at unitarity. The inset shows the red shift of the 
total dynamic structure factor with decreasing temperature. 

mentally, Sfj,(q, ui) can be measured by a proper choice 
of the detuning of the laser beams that results in different 
couplings to the two spin components [4]. In Fig. 2, we 
show the evolution of S a a' (q, u) with increasing interac- 
tion strength (from bottom to top). At lower tempera- 
tures (left column), ^^(q, w) grows and becomes com- 
parable to Sff(q,u;) at 1/kpa = 1.0. It should be noted 
that our results at T = 0.5Tp are qualitatively reliable 
and are presented for illustrative purposes only. We ex- 
pect the predictions at T = Tp to be more quantitative, 
as estimated conservatively from the virial expansion of 
the equation of states for a trapped strongly interacting 
Fermi gas [12]. 

Fig. 3 gives the temperature dependence of the DSF at 
unitarity. As expected, S'-f^q, uS) increases rapidly with 
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Figure 4: (Color online). Dependence of static structure fac- 
tors on the momentum (a) and interaction strength (b) at 
T = T F . 

decreasing temperature. As a consequence, the peak of 
total dynamic structure factor shifts towards the molec- 
ular recoil frequency, as indicated in the inset. This red- 
shift of atomic peak was indeed observed experimentally 
for a unitarity Fermi gas at T ~ OATp [13 ]. 

We have so far restricted ourselves to cases with large 
q. The momentum dependence can be conveniently 
illustrated by the static structure factor S aa i (q) = 
(2/N) J du)S aa >(q,uj). In Fig. 4, we present .Sw(q) 
as a function of momentum (Fig. 4a) and interac- 
tion strength (Fig. 4b) at the Fermi temperature Tp. 
The spin-antiparallel structure factor depends strongly 
on momentum and interaction strength. In contrast, 
the spin-parallel structure is nearly always unity, partly 
due to the autocorrelations among identical spins. To 
better understand this, we consider the short-range be- 
haivor of the spin-parallel structure factor in real space, 
S aa (r, r') = (2/JV) {h a (r) n a (r')}. For parallel spins, 
there is no singularity as r' — > r. Thus, to a good approx- 
imation S aa ( r , r ') — (2/A) {ha (r)) 8 (r — r'), where the 
delta function comes from the anticommutator of fermion 
operators. After the Fourier transformation and average 
over the trap, we obtain S^^q) — 1. 

CONCLUSIONS 

In conclusion, we have developed a quantum virial ex- 
pansion for dynamical properties of strongly correlated 
systems and have computed the dynamic density re- 
sponse of a strongly interacting Fermi gas at tempera- 
tures T ~ Tp. The experimentally observed transition 
from atomic to molecular response at low temperatures 
(~ Q.lTp) has been reproduced at a qualitative level. We 
anticipate that future Bragg spectroscopy at high tern- 
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peratures will lead to a quantitative agreement. Alter- 
natively, by including higher-order virial expansion func- 
tions, we may extend the validity of our results closer 
to the characteristic critical temperature ~ 0.2Tp [l^ . 
We emphasize that our virial theory is efficient for in- 
vestigating other basic dynamical properties, such as the 
spectral function of single-particle Green function. In 
this respect, it may shed light on solving the paradox of 
the pseudogap phenomenon at unitarity (l5j . We leave 
this possibility in a future publication. 
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